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The  report  describes  the  results  of  Vtwo  separate  experimental  investigations. 
The  first  concerns  the  electrical  excitation  of  the  HC1  chemical  laser  and  was 
motivated  by  the  potential  improvement  of  the -laser  performance  by  the  vibrational 
enhancement  of  the  reaction  Cl  + H2(v=1 ) s*  HC1 (v=v,l ) + H.  Detailed  kinetic 
modelling  of  the  e-beam  and  e-beam  sustained  discharge  experiments  was  used  in 
the  interpretation  of  the  experimental  data.  The  rate  constant  for  the  vibra- 
tional ly  enhanced  reaction  was  40  times  larger  than  for  the  reaction  with 
tf£(v=0)  at  T - 30Q/K.  This  enhanced  rate  is  not  fast  enough  to  lead  to 
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20.  Abstract  (continued) 


efficient  use  of  the  Hp( v=l ) excited  by  electron  impact.  The  maximum  electrical 
to  laser  conversion  efficiency  obtained  was  3 percent. 

The  second  part  of  the  report  describes  measurements  of  the  threshold  pump- 
ing power  density  in  e-beam  sustained  discharge  excitation  of  KrF  and  XeF  lasers. 
Stable  discharges  were  obtained  in  Ar-Kr-F^and  Ar-Xe-NFj  mixtures' at  1 atm  for 
times  greater  than  0.5  psec  and  lasing  pulse  width  as  long  as  0.5  mSoc.  The 
threshold  pump  power  densities  varied  inversely  with  the  cavity  build  up  time. 
Based  on  the  present  results  it  is  evident  that  efficient  laser  operation  in 
KrF  and  XeF  could  be  achieved  in  a 1-m  long  sustained  discharge  with  input  power 
densities  in  the  range  of  50  to  100  kW/cm3  over  a 1 psec  time  scale. 
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PREFACE 


This  report  is  presented  in  two  parts.  The  first  concerns 
the  investigation  of  e-beam  controlled  discharge  initiation  of  the 
HC1  chemical  laser  system.  The  second  part  describes  an  experi- 
mental investigation  of  the  threshold  power  density  requirements 
for  e-beam  sustained  discharge  excitation  of  XeF  and  KrF. 

This  research  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by  ONR  under 
Contract  No.  N00014-75-C-0052. 
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SECTION  I 


INTRODUCTION 

The  rate  of  an  endothermic  chemical  reaction  can  be  significantly 

1 2 

increased  by  vibrational  excitation  of  the  reactant  molecules.  * For  the 

near  thermoneutral  reaction  Cl  + H2  -*■  HC1  (AH  * 1 kcal),  theoretical  and 

experimental  estimates  have  been  made  for  the  accelerated  reaction  rate  of 

3 4 

vibrationally  excited  H~.  Estimates  range  from  a factor  of  25  ’ to 

5 c 

2000  for  the  ratio  of  rate  constants  ku  (v=l)/ku  (v=0).  If  the  larger 

T12  Ho 

factor  were  correct,  substantial  improvements  could  be  expected  in  the  per- 
formance of  the  electrically  excited  HC1  chemical  laser.  The  reason  for 
this  improvement  can  be  understood  by  considering  the  rate  data  presented 
in  Table  1.  Reactions  (1)  and  (2)  constitute  the  HC1  chain  reaction  steps 
leading  to  an  HC1  laser.  Reaction  (2)  is  the  slower  rate-controlling  step 
which  must  regenerate  H atoms  for  the  highly  exothermic  "hot"  reaction  (1). 
Laser  operation,  therefore,  dictates  a high  initial  H or  Cl  atom  density 
which  could  produce  rapid  V-T  deactivation  of  excited  HC1  molecules  and  re- 
duced laser  efficiency.  If,  on  the  other  hand,  reaction  (3)  were  three 
orders  of  magnitude  faster  than  (2),  100  times  smaller  concentrations 
of  Cl  atom  would  be  required  at  modest  hydrogen  vibrational  temperatures 
(2000  °K)  to  produce  the  same  total  rate  of  H-atom  production  as  reaction 
(2). 


The  goal  of  the  present  study  was  to  investigate  the  performance 
of  an  HC1  chemical  laser  making  use  of  reaction  (3)  and  to  evaluate  the 
magnitude  of  the  rate  enhancement.  A short  pulse  (1  usee)  electron  beam 
stabilized  discharge  has  been  used  to  initiate  the  HC1  chemical  laser  in 
mixtures  of  Ar,  H2,  and  Cl2.  Reliable  estimates  could  be  made  for  the  degree 
of  H2  vibrational  excitation  since  the  electron  energy  distribution  is  fairly 
well  characterized  for  these  externally  ionized  discharges.  A comprehensive 
computer  model  of  the  HC1  laser  and  chemical  kinetics  used  the  estimates 
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Table  1 

Selected  Reaction  Rate  Constants  for  HC1  Chemical  Laser  (T  = 300  °K) 


Reaction 


1.  H + Cl2  -^HCl  (v=l  ,2,3, ... ) + Cl 

l< 

2.  Cl  + H2  -co1dr  HC1  (v=0)‘+  H 

3.  Cl  + H2  (v=1  ) ->  HC1  (v'=0,l ) + H 

4.  Cl  + HC1  (v=l)  HC1  (v  + 0)  + Cl 

5.  H + HC1  (v=l)  - HC1  (v=0)  + H 


Rate  Constant 
(cm3/mole-sec) 

3.8  x 1012 
1.5  x 1010 
6.3  x 1011 

4.8  x 1012 
4.0  x 1012 
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for  hydrogen  vibrational  excitation  to  calculate  the  energies  and  intensi- 
ties of  the  various  HC1  vibrational  laser  transitions  and  time-dependent 
concentrations  of  the  excited  HC1  levels  (v=0, 1,2,3),  H and  Cl  atoms.  Cl 2 
and  H2  (v=0,1). 

Three  principal  diagnostic  tools  were  used  to  study  the  HC1  chemical 

laser  under  a variety  of  discharge  conditions.  Time-resolved  spectroscopic 

measurements  were  made  of  the  overall  pulse  shape  and  the  individual  HC1 

35  37 

laser  lines,  for  both  H Cl  and  H Cl  isotopes.  Consumption  of  the  Cl2  was 
monitored  using  visible  absorption  spectroscopy.  Laser  output  energy  mea- 
surements were  also  made  as  an  additional  constraint  for  evaluating  the 
vibrationally  enhanced  reaction  rate  constant  with  the  computer  model. 

Good  agreement  between  the  experimental  measurements  and  the  computer 
simulation  has  been  achieved  and  will  be  described  more  fully  in  the  main 
text.  A value  of  40  for  the  vibrational  rate  enhancement  (k^(v=l )/k^(v=0)) 
has  been  found  to  be  consistent  with  the  experimental  observations. 

Sections  II  and  III  describe  the  experimental  apparatus,  the  oper- 
ating characteristics  of  the  HC1  chemical  laser,  and  the  diagnostic  measure- 
ments. A computer  model  of  the  reaction  kinetics  and  laser  characteristics 
are  presented  in  Section  IV  together  with  comparisons  with  experimental 
measurements. 


SECTION  II 


DESCRIPTION  OF  EXPERIMENT 

Figure  1 is  a general  schematic  diagram  of  the  laser  excitation 
apparatus  and  diagnostics.  The  laser  cavity  was  formed  by  two  2-m  radii - 
of-curvature  totally  reflective  mirrors  (R).  CaF^  flats  were  used  as 
laser  windows  (W)  and  output  coupler  (0).  Infrared  detector  (I)  and  mono- 
chromator (M)  measured  the  time-resolved  spectroscopic  content  of  the 
laser  pulse.  An  energy  meter  (E)  and  infrared  detector  (I')  were  used 
at  other  times  to  measure  the  total  laser  energy  and  pulse  shape.  A short- 
pulse  (1  ysec)  cold-cathode  electron  beam  source  was  used  to  stabilize 
discharges  in  mixtures  of  Ar,  F^,  and  C^.  A cross-sectional  drawing  of 
the  essential  features  of  the  discharge  chamber  and  electron-beam  emitter 
structure  is  shown  in  Figure  2.  The  electron  gun  operated  at  160  kV  and 
delivered  3.5  A/cm2  through  the  1-mil  aluminized  Kapton  foil  which  separated 
the  discharge  chamber  from  the  diode  vacuum  chamber.  By  adjusting  the 
spacing  from  the  graphite  (or  thin  tantalum  foil)  emitter  to  the  electron- 
beam  anode  screen,  the  current  pulse  amplitude  and  duration  could  be  altered. 
Experiments  were  successfully  run  with  10-usec  pulse  durations  at  0.3  to 
0.4  A/cm2  current  densities. 

The  discharge  chamber  consisted  of  a high-density  polyethylene  body 
which  provided  high-voltage  electrical  isolation  for  the  aluminum  anode  as 

shown  in  Figure  2.  An  array  of  10  coaxial  cables  connected  the  discharge 
anode  to  a triggered  gas  switch  and  2-pF  energy  storage  capacitor.  Elec- 
tronic delay  circuits  allowed  the  discharge  to  be  fired  at  any  time  with 
respect  to  a trigger  pulse  which  initiated  the  electron  gun.  The  discharge 
current  rise  time  was  typically  200  to  300  nsec  and  was  limited  by  the  in- 
ductance of  the  connecting  cables  and  capacitor  (approximately  100  nH). 

The  optical  aperture  of  the  discharge  region  is  reduced  by  a 2-cm  diameter 
channel  in  the  high  density  polyethylene  at  each  Brewster  angle  window  as 
shown  in  Figure  1.  Figure  2 schematically  shows  the  effective  system  optical 
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Figure  1.  Schematic  Diagram  of  the  HC1  Laser  Excitation  Apparatus 
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aperture  with  respect  to  the  5-cm  wide  discharge  anode.  The  discharge  anode 
was  50  cm  long,  and  the  overall  length  of  the  chamber  was  75  cm.  Chlorine 
absorption  measurements,  to  be  described  later,  indicated  very  little  C 1 ^ 
disappearance  for  5-or  10-minute  sample  storage  times  in  the  cavity.  The 
white  polyethylene  walls  adjacent  to  the  discharge  region  did  discolor  after 
several  hundred  shots,  partly  due  to  radiation  damage  from  the  electron  beam. 

The  argon,  hydrogen,  and  chlorine  gases  were  obtained  from  Matheson 
and  were  either  pre-puri fied  grade  (Ar,  H^)  or  high  purity  grade  (Cl^).  and 
were  used  without  further  purification.  The  stated  purities  were  99.998%  for 
Ar,  99.95%  for  ^ and  99.5%  for  C^.  Matheson  flowmeters  were  used  to  adjust 
the  relative  concentrations  of  the  gases.  The  meters  were  connected  together 
and  the  gases  flowed  through  13  ft  of  copper  tubing  (1/4-in.  0D,  1/8-in.  ID) 
for  good  mixing  before  entering  the  discharge  chamber.  There  was  a 20  to  30 
percent  difference  in  the  laser  output  energies  observed  under  static  and 
flowing  gas  conditions.  The  overall  laser  spectrum  and  pulse  shape  were 
essentially  the  same  for  both  conditions.  All  the  chlorine  visible  absorp- 
tion data  were  taken  using  flowing  gas  mixtures. 

Spectroscopic  measurements  of  the  individual  HC1  laser  transitions 

were  made  with  quarter-m»ter  (Jarrel-Ash)  and  one-meter  (McPherson)  grating 

monochromators.  The  frequency  resolution  of  the  smaller  instrument  was 

± 0.8  cm-1.  Higher  resolution,  ± 0.2  cm-1,  was  obtained  using  the  one-meter 

37 

monochromator,  thereby  permitting  unambiguous  identification  of  the  H Cl 

35  37 

lines.  Individual  H Cl  and  H Cl  P-branch  transitions  were  recorded  in  a 
shot-by-shot  search  of  their  respective  vibrational  bands.  A cooled  (77  °K) 
gold-doped  germanium  detector  (Raytheon,  QKN  1568),  100-nsec  rise 
time,  was  used  for  all  the  infrared  laser  measurements.  Matching  100  fi 
coaxial  cable  and  resistive  terminator  were  used  for  optimum  time  response. 

Figure  3 shows  the  experimental  optical  arrangement  used  for  the 
chlorine  absorption  measurements.  A high  brightness,  100-watt  mercury  arc 
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lamp  (Osram,  HBO  100  W/2)  provided  an  intense  365  nm  source.  Four  Corning 
glass  filters  (two  each  5860  and  5840)  were  used  to  isolate  the  365-nm 
spectral  region  emitted  by  the  source.  Care  was  taken  to  further  attenuate 
the  light  emitted  by  the  mercury  arc  using  screens  so  that  significant 
chlorine  dissociation  did  not  occur.  Two  criteria  were  established  for  de- 
termining the  maximum  light  flux  which  could  be  transmitted  through  the  gas 
samples.  First,  the  measured  absorption  did  not  decrease  over  a 10-  or  15- 
minute  time  span  with  the  gas  sample  slowly  flowing  through  the  discharge 
chamber.  As  a second  check  for  pre-reaction,  the  chlorine  mole  fraction 
measured  by  the  visible  absorption  technique  should  agree  within  10  percent 
with  the  concentration  found  using  the  relative  filling  time  of  the  dis- 
charge chamber  (typically  over  the  range  100  to  300  torr)  with  each  of  the 
three  component  gases.  Since  the  mercury  arc  operated  DC  and  the  experimental 
measurements  were  pulsed,  a chopper  wheel  with  a single  1 in.  opening, 
operating  at  1725  rpm,  provided  additional  attenuation  by  reducing  the  duty 
cycle  of  the  arc  illumination  of  the  sample.  A photodiode  viewing  the 
light  reflected  from  the  glass  filters  provided  a synchronization  signal  for 
the  electron  beam  and  discharge  circuitry.  A zero  percent  baseline  absorp- 
tion was  measured  for  each  experiment  by  shutting  off  the  chlorine  flow 
following  an  excitation  pulse,  waiting  3 to  4 minutes  for  the  flowing  argon 
and  hydrogen  to  purge  the  chamber  of  chlorine,  and  then  recording  the  light 
signal  on  the  oscilloscope.  One-hundred  percent  absorption  was  also  re- 
corded on  the  same  oscilloscope  trace  by  physically  blocking  the  light  reach- 
ing the  detection  system. 

Two  problems  encountered  in  these  measurements  were  the  light  emitted 
by  the  discharge  itself  and  the  electrical  noise  generated  by  the  electron 
beam  and  discharge.  The  latter  difficulty  was  remedied  using  a shielded 
and  electrically-isolated  Faraday  cage  for  the  photomultiplier  and  its 
associated  electronics  (power  supply  and  oscilloscope).  A quarter-meter 
monochromator  was  used  with  the  photomultiplier  to  alleviate  the  background 
luminosity  problem.  Attenuation  of  the  mercury  arc  source  was  also  adjusted 
for  maximum  light  flux  to  overcome  background  emission  from  the  discharge 
chamber. 
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A Gen-Tec  energy  meter  (ED-200)  was  used  to  measure  the  total  single- 
pulse  HC1  laser  energy.  When  the  energy  meter  was  used,  an  infrared  detec- 
tor monitored  the  laser  pulse  time  history  (Fig.  1).  A tilted  CaF^  flat 
inside  the  laser  cavity  provided  approximately  8 percent  output  coupling 
(round-trip)  at  30  degrees  incidence  angle.  The  total  cavity  losses  were 
estimated  to  be  13  percent  when  mirror  reflectivity  losses  and  scattering 
from  window  surfaces  were  included. 


SECTION  III 


EXPERIMENTAL  RESULTS 


A.  Laser  Pulse  Shape  and  Spectroscopic  Analysis 

HC1  laser  pulse  time  histories  for  three  different  excitation  condi- 
tions are  shown  in  Figure  4.  The  gas  composition  (Ar/H2/Cl2  49.1/49.1/1.8) 
and  pressure  (400  torr)  were  the  same  for  all  three  examples  presented. 

The  electron  beam  and  electric  discharge  initiation  pulses  began  at  t = 0 
and  were  terminated  at  t = 1 usee.  Using  the  electron  beam  alone  for 
initiation  of  the  chemical  laser  always  produces  a pulse  with  three  promi- 
nent peaks  as  seen  in  Figure  4(a).  Two  different  values  of  E/N  were  ex- 
tensively studied.  A low  value  (1.55  x 10~16V-cmz)  was  found  which 
significantly  enhanced  the  laser  output  compared  to  the  electron  beam 
alone  for  initiation.  The  high  E/N  value  (2.5  x 10'16)  was  the  maximum 
voltage  the  appratus  could  reliably  withstand  without  arc  formation  in 
the  laser  gas.  The  pulse  shape  of  the  low  E/N  discharge-initiated  HC1 
laser  was  also  quite  characteristic  and  reproducible.  High  E/N  discharges 
produced  laser  pulses  which  did  vary  somewhat  from  day-to-day,  although 
the  overall  characteristics  were  a short  duration  (typically  2.1  usee)  and 
smooth,  steep  drop-off  after  the  laser  reached  maximum  intensity.  On  some 
traces,  three  small  distinct  peaks  spaced  by  approximately  150  nsec,  or 
two  peaks,  spaced  by  300  nsec,  were  seen  superimposed  on  the  general  peak 
intensity  curve,  for  the  high  E/N  laser  pulses.  Some  of  this  variability 
in  the  appearance  of  the  peak  laser  intensity  was  correlated  with  slight 
variations  (100  to  200  nsec)  in  the  electron  beam  pulse  duration,  which 
is  essentially  the  same  as  the  discharge  pulse  duration. 

Longest  duration  lasing  was  observed  for  low  E/N  or  low-energy  input 
discharges  at  400  torr.  The  slightly  earlier  laser  threshold  time  of  the 
high  E/N  case  (Figure  4(c))  can  be  qualitatively  explained  by  the  higher 
H-atom  and  H2(V)  production  rates  in  the  more  energetic  discharge.  Molecular 
chlorine  consumption  by  the  chemical  reaction  and  buildup  of  chlorine 
atoms  causes  the  pulse  to  end  quite  abruptly.  The  enhanced  cold  reaction 
rate  caused  by  hydrogen  vibrational  excitation  is  qualitatively  evident  in 
Figure  4(b)  as  the  low-intensity,  long  duration  tail  of  the  laser  pulse. 


RELATIVE  INTENSITY  (ARB.  UNITS) 
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Figure  4.  Experimentally  Measured  HC1  Laser  Pulse  Shapes  for  Electron 
Beam  and  Discharge  Excitation  (1-psec  Pulse  Duration).  Gas 
pressure  and  mixture  were  400  torr  and  Ar/H2/Cl2  49/49/2. 
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Detailed  spectroscopic  measurements  were  made  for  the  three  cases 
of  Figure  4,  and  this  information  is  graphically  presented  in  Figures  5,  6, 
and  7.  Laser  emission  from  3-2  transitions  was  generally  much  weaker  than 
2-1  and  1-0  bands.  Only  in  the  electron  beam  initiation  case  (Fig.  5)  did 
the  1-0  transitions  reach  threshold  significantly  later  than  the  3-2  and 
2-1  bands.  The  late-time  emission  in  the  low  E/N  discharge  experiments 
was  dominated  by  1-0  transitions. 

37 

Individual  H Cl  laser  lines  were  identified  and  time  resolved  using 

35 

the  one-meter  monochromator.  The  average  spacing  of  a pair  of  H Cl  and 
37 

H Cl  lines  is  1.8  cm-1,  which  is  essentially  the  resolution  limit  of  the 
smaller  quarter-meter  instrument.  However,  accidental  near  coincidences, 
such  as  2-1  P ( 7 ) H37C1  and  1-0  P(ll)  H35C1  (Av  * 0.6  cm*1),  were  only  re- 
solvable using  the  larger  instrument.  Under  low  resolution,  some  of  the 
P(J)  laser  pulses  in  Figures  5,  6,  and  7,  for  example  1-0  P(7)  (Fig.  6), 

would  appear  to  have  two  time-resolved  peaks,  whereas  in  fact  the  second 

37 

pulse  was  the  time-delayed  H Cl  peak. 

37 

Lasing  from  H Cl  2-1  and  1-0  transitions  affects  the  laser  pulse 
shape  as  can  be  seen  in  Figure  8,  where  a composite  laser  pulse  has  been 
constructed  from  the  data  of  Figure  6.  The  agreement  with  the  experimental 
pulse  shape  (Figure  4(b))  is  quite  good. 

35  37 

As  a result  of  the  relative  isotopic  abundance  of  Cl  and  Cl,  the 

Inherent  gain  of  an  H37C1  line  is  one-third  that  of  the  same  H35C1 

37 

vibrational-rotational  transition.  Therefore,  certain  H Cl  transitions 

35 

were  completely  missing,  even  though  the  corresponding  H Cl  transitions 

37 

were  observed.  The  observed  H Cl  transitions  had  a longer  time-to-threshold 
and  approximately  one-third  lower  intensity  than  the  identical  transitions 
in  H35C1. 

37 

Since  the  H Cl  transitions  had  lower  gain  and  required  longer  time 
to  reach  threshold,  the  upper  vibrational  levels  were  much  more  sensitive 
to  the  chemical  pumping  reactions  and  vibrational  decay  process.  Therefore, 
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Figure  5.  Time  Histories  and  Relative  Intensities  of  the  HC1  Laser 

Transitions  Observed  Using  Only  the  Electron  Beam  for  Exci- 
tion.  Gas  pressure  was  400  torr  and  the  gas  mixture  was 
Ar/H2/Cl2  (49/49/2). 
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Figure  8.  Reconstructed  HC1  Laser  Pulse  Shape  Using  the  Low  E/N  Discharge 
Data  Presented  in  Figure  6 
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37  35 

the  information  obtained  from  the  time-resolved  H Cl  and  H Cl  spectro- 
scopic measurements  provided  a sensitive  internal  standard  and  cross-check 
for  developing  the  computer  model  described  in  Section  IV. 

B.  Chlorine  Visible  Absorption  Measurements 

Measurements  of  the  time-dependent  molecular  chlorine  concentration 
were  made  for  various  excitation  conditions.  Figure  9(a)  shows  an  absorp- 
tion trace  for  a high  E/fl  discharge.  The  early  spike  near  1 psec  was  caused 
by  light  emission  from  the  discharge.  Estimates  for  the  direct  dissociation 
of  the  chlorine  by  the  electron  beam  and  electrical  discharge  could  be  made 
from  the  initial  drop  in  365  nm  absorption  iimediately  following  the  exci- 
tation pulse.  However,  significant  chlorine  dissociation  only  followed  a 
high  E/N  discharge.  Therefore,  the  most  meaningful  use  of  the  initial 
chlorine  absorption  data  was  to  establish  an  upper  bound  to  the  direct 
chlorine  dissociation.  Plotting  the  chlorine  concentration  versus  time  for 
the  linear  portion  of  the  absorption  decay,  as  in  Figure  9(b),  yielded  a 
global  value  for  the  chlorine  reaction  rate. 

The  sensitivity  of  the  chlorine  absorption  measurements  was  reduced 
by  the  presence  of  a 12.75  cm  region  of  unexcited  gas  at  each  end  of  the 
laser  cavity.  The  total  optical  path  for  the  absorption  measurements  was 
75.5  cm,  with  25.5  cm  containing  unexcited  gas.  The  absorption  coefficient 
of  Cl ^ at  365  nm  is  3.73  x 10‘3  cm-1  torr"1.®  Hence,  for  a typical  1.75 
percent  chlorine  mixture  at  400  torr,  50  percent  absorption  will  always  be 
present  in  the  window  regions.  This  effect  is  evident  in  Figure  9(a) 
where  the  absorption  asymtotical ly  approaches  50  percent.  After  several 
milliseconds,  the  chlorine  absorption  in  the  end  regions  also  decreases, 
corresponding  to  the  propagation  of  the  reaction  into  these  regions. 

C.  Laser  Energy  Measurements 

Parametric  laser  energy  measurements  were  made  for  a variety  of  gas 
mixtures  and  discharge  conditions.  Figure  10  shows  the  laser  energy  as  a 
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Figure  TO.  Measured  HC1  Laser  Energy  as  a Function  of  Cl p and  H~  Mole 
Fractions  Using  Electron  Beam  Excitation 
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function  of  hydrogen  and  chlorine  concentrations  using  only  the  electron 
beam  to  initiate  the  chemical  laser.  The  laser  output  was  fairly  insensi- 
tive to  the  chlorine  and  hydrogen  concentrations  over  the  range  shown  in 
Figure  10.  However,  Figure  11  shows  that  for  a constant  discharge  energy 
input  and  E/N,  there  is  a strong  dependence  of  the  laser  energy  on  the 
chlorine  fraction.  The  energy  input  to  the  discharge  was  kept  constant  by 
attenuating  the  electron  beam  current,  hence  maintaining  constant  discharge 
current,  as  the  chlorine  fraction  was  reduced.  This  adjustment  in  the  elec- 
tron beam  current  is  necessary  because  of  the  large  (-1 0~ 1 °cm3/sec)  electron 
attachment  rate  of  the  chlorine  molecule.  The  laser  energy  using  only  the 
electron  beam  for  excitation  is  also  shown  in  Figure  11.  In  contrast  to 
Figure  10,  the  electron  beam  current  is  not  the  same  for  each  chlorine  frac- 
tion studied. 

The  highest  electrical  efficiency  obtained  was  4.5  percent  in  a 
high  Cl 2 mole  fraction  (3.5  percent)  mixture  as  shown  in  Figure  11. 

Table  2 summarizes  the  laser  output  energy  and  electrical  efficiency 
for  the  excitation  conditions  shown  in  Figure  4.  The  laser  energy  was 
approximately  ten  times  higher  using  the  discharge  than  it  was  using  the 
electron  beam  alone.  Higher  E/N  discharges  would  produce  more  H atoms  (see 
Sec.  IV)  which  would  improve  the  laser  output  by  accelerating  the  rate  of 
reaction  (1)  (Table  1)  but  reduce  the  electrical  efficiency.  The  quantum 
efficiency  of  producing  an  H atom  is  at  best  10  percent,  based  on  a 9-eV 
H^  dissociation  energy  and  assuming  production  of  HC1  (v=2),  and  ignoring 
V-T  deactivation. 

As  expected  for  a laser  system  which  is  strongly  dominated  by  V-T 
loss  mechanisms,  the  input  energy  pulse  duration  plays  a crucial  role.  The 
cold  cathode  electron  gun  and  discharge  could  be  operated  at  excitation 
pulse  lengths  as  long  as  15  usee  or  as  short  as  several  hundred  nsec.  In 
the  course  of  early  studies,  it  was  found  that  by  decreasing  the  pulse 
duration  from  10  to  5 usee,  for  the  same  total  electrical  energy  input, 
the  laser  energy  doubled.  These  findings  led  to  the  intensive  studies  con- 
ducted with  the  1-usec  excitation  pulse  configuration.  For  pulse  durations 
less  than  1 usee,  the  laser  energy  scaled  directly  with  the  energy  input 
pulse  length. 
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Figure  11.  Measured  HC1  Laser  Energy  as  a Function  of  Cln  Mole  Fraction 
Using  Discharge  Excitation  at  Constant  Energy^Input. 
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Table  2 

Summary  of  Laser  Output  Energies  and  Electrical  Efficiencies 


Excitation1  ^neroy 

Conditions  (mJ) 

3 

Electrical 
Energy 
Input  (J) 

Electrical 

Efficiency 

(*) 

Chemical4 

Efficiency 

(%) 

5 

Electron  Beam 
(1  ysec  pulse) 

3.5  A/cnr  20  - 30. 

1.2 

~ 2.0 

~ 0.2 

Electron  Beam + Low 
(1 .6  x 10-16V-cm2) 

E/N  Discharge 

(1  ysec  pulse)  200  - 230 

8.1 

~ 2.6 

~ 1.9 

Electron  Beam + High 
(2. 5 x 10-1 6V-cnr ) 

E/N  Discharge 

(1  ysec  pulse)  244  - 260 

12.6 

~ 2 

~ 2.2 

^Gas  mixture  Ar/Hg/C^  (49/49/2)  at  400  torr. 

2 

A range  of  laser  energies  is  given  which  covers  identical  dis- 
charge conditions,  data  taken  on  different  days. 

3 

162  cm3  active  volume  was  viewed  by  laser  optics. 

4 

Based  on  44  kcal/mole  of  Cl 2 present. 

5 

The  Ion  pair  production  rate  corresponding  to  3.5  A/cm?  delivered 
by  the  electron  gun  was  1.6  x 10?1  sec-1  and  28  eV  per  ion  pair  was 
deposited  in  the  gas. 


SECTION  IV 

COMPUTER  MODEL  DESCRIPTION  AND  COMPARISONS  WITH  EXPERIMENT 


The  computer  model  was  divided  into  two  independent  calculations. 
Purely  discharge  and  electron  impact-related  phenomena,  such  as  vibrational 
excitation  of  hydrogen,  were  treated  in  a Boltzmann  calculation  of  the 
electron  distribution  function.  The  second  set  of  calculations  dealt  with 
the  reaction  kinetics  and  laser  gain  and  used  the  fractional  distribution 
of  electron  beam  and  discharge  energy  into  various  molecular  excitation 
channels  as  input  data.  Subsections  A and  B describe  the  essentials  of  the 
electron  excitation  calculations  and  the  overall  features  of  the  chemical 
kinetics  and  laser  gain  program.  A set  of  chemical  rate  constants  and 
values  for  the  fractional  distribution  of  electrical  energy  input  were  found 
which  provided  good  agreement  with  the  experimentally  measured  laser  and 
kinetic  properties.  These  parameters  will  be  discussed  in  Subsection  C. 


A.  Analytical  Model  of  Electron  Beam  and  Electric  Discharge  Physics 


e-Beam  Physics 


The  production  rate,  S,  of  secondary  electrons  from  the  high-energy 
primary  electrons  (160  keV), 


ep  + Ar  ->  Ar+  + ep  + es 

can  be  determined  from  the  e-beam  current  density,  mass  density,  effective 
ionization  energy,  and  primary  electron  energy  loss  rate  for  each  species  of 
the  gas  mixture.  Effective  ionization  energies  for  and  Ar  were  obtained 
from  Whyte,^  and  the  electron  energy  loss  rates  were  obtained  from  Berger 
and  Seltzer.®  It  has  recently  been  demonstrated  by  Searles  and  Hart^  ^ 
that  the  energy  deposition  rates  calculated  using  the  Berger  and  Seltzer 
tables  are  two  to  four  times  smaller  than  that  obtained  from  calculations 


25 


which  consider  scattering  of  the  incident  high  energy  electrons  by  the  foil 
and  multiple  scattering  in  the  gas.  Hence,  the  ion-pair  production  rate  S 
calculated  using  the  Berger  and  Seltzer  tables  (8  x 1020/cm3-sec)  for  the 
standard  gas  mixture  studied,  49%  Ar,  49%  H2,  2%  Cl^,  was  multiplied  by 
2 to  yield  1.6  x 1021/cm3-sec.  This  value  of  the  ion-pair  production  rate 
was  found  to  be  consistent  with  the  observed  discharge  current  and  atom 
(H  and  Cl)  production  rates,  as  will  be  discussed  below. 

Because  of  the  large  electron  attachment  rate  coefficient  for  Cl 2 » 
the  total  attachment  rate  can  be  equated  to  the  ion-pair  production  rate 
and  dissociative  recombination  processes  can  be  ignored.  Hence,  the  ob- 
served discharge  current  only  depends  upon  the  electron  drift  velocity  and 

attachment  rate  coefficient.  For  Cl0»  the  published  thermal  electron  attach- 

^ (12)  (13) 

ment  rate  constants  range  from  3 x 10_loV  ' to  3.7  x 10"9V  'cm3/molecule- 
sec.  Very  little  information^  exists  on  the  electron  energy  dependence 
(above  0.25  eV)  of  the  rate  coefficient.  From  the  present  measurements  of 
the  discharge  current  density  (24  A/cm2)  in  low  E/N  (1.6  x 10" 16  V-cm2) 
discharges,  in  Ar/I^/C^  (49/49/2)  mixtures,  an  electron  attachment  rate 
coefficient  of  1.5  x 10“10  cm3/molecule-sec  for  C 1 2 has  been  found  using  a 
calculated  drift  velocity  of  3.3  x 106  cm/sec.  The  procedure  for  calculating 
the  drift  velocity  will  be  outlined  in  the  discharge  physics  section  below. 

Some  elementary  time-dependent  computer  calculations  were  made  for 

15 

the  electron  detachment  process  H + Cl  -*■  HC1  + e,  using  the  published 
rate  constant  (9  x 10"10cm3/molecule-sec)  for  this  reaction.  In  the 
absence  of  the  detachment  reaction,  the  electron  density  reached  the 
steady  state  value  (3  x 1013/cm3)early  in  the  pulse  (100  nsec),  but  then 
decayed  to  half  this  value  by  the  end  of  the  1 psec  pulse.  The  inclusion 
of  the  H-atom  detachment  reaction  prevented  the  decay  of  the  electron 
density  over  the  entire  pulse  duration.  However,  the  H atom  reaction 
did  not  significantly  affect  the  overall  magnitude  of  the  electron 
density,  which  was  still  strongly  determined  by  the  Cl 2 attachment  rate 
coefficient. 
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Discharge  Physics 

In  the  presence  of  an  electric  field,  it  is  necessary  to  solve  the 
Boltzmann  equation  for  the  electron  distribution  function  in  order  to  obtain 
the  drift  velocity,  mean  electron  temperature,  and  excitation  rates  for  the 
various  electron  impact  inelastic  collision  processes  involved. 

For  a given  gas  mixture  and  applied  electric  field  normalized  by  the 
gas  density  (E/N),  the  Boltzmann  transport  equation  can  be  solved  if  the 
elastic  and  inelastic  electron-impact  cross  sections  are  known  as  a func- 
tion of  elect/on  energy.  These  cross  sections  are  well  known  for  the 
dominant  gases,  and  Ar,  in  the  gas  mixture  under  study  but  are  not  avail- 
able for  Cl 2 (except  for  the  dissociative  attachment  process). 

Cross  sections  for  argon  were  obtained  from  Engelhardt  and  Phelps. ^ 

17 

For  hydrogen,  the  momentum  transfer  cross  section  came  from  Gibson,  the 

18 

rotational  and  vibrational  excitation  cross  section  from  Linder  and  Schmidt, 
and  the  dissociation  cross  section  from  Chun,  Lin  and  Lee.^  The  drift 
velocity  and  characteristic  electron  energy  are  shown  as  a function  of  E/N 
in  Figure  12  for  a 50%  H2-50%  Ar  mixture.  Fractional  electrical  power  dis- 
tribution as  a function  of  E/M  is  shown  in  Figure  13. 

A check  on  the  self-consistency  of  the  model  for  both  e-beam  and  dis- 
charge properties  was  made  by  comparing  the  experimentally  measured  discharge 
current  with  that  given  theoretically.  The  electron  density  was  calculated 
by  equating  the  ion-pair  production  rate  and  the  Cl2-dominated  electron 
attachment  rate  (S  = 6 N&  Nci^*  The  electron  drift  velocity  was  obtained 
from  the  Boltzmann  code  calculations.  The  discharge  current  was  then  found 
by  the  product  of  electron  density,  drift  velocity  and  electron  charge.  The 
agreement  was  generally  within  10  to  20  percent,  which  is  quite  satisfactory. 

Since  data  for  the  electron  impact  dissociation  cross  section  of  Cl2 
was  not  available,  parametric  values  for  the  fraction  of  discharge  energy 
entering  Cl ^ dissociation  were  tried.  It  was  found  that  only  10  percent 
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or  less  of  the  discharge  energy  needed  to  be  placed  in  C 1 2 dissociation  to 
give  good  agreement  with  the  observed  C 1 2 burn  rate  and  laser  properties 
discussed  in  Section  II.  This  implies  4.5  percent  dissociation  of  the  Cl^ 

by  the  discharge  (E/N=  1.6x  1 0~ 1 6 V-cm2).  For  comparison  purposes,  1%  F? 

20  L 
dissociation  has  been  measured  for  an  electron  beam-stabilized  discharge 

in  mixtures  of  He/F^^.  Reference  20  also  proposes  a mechanism  involving 
heating  of  F~  ions,  by  the  discharge  field,  to  temperatures  of  the  order 
of  1 eV.  At  this  energy,  the  electron  can  rapidly  detach  from  F"  in  heavy 
particle  collisions,  and  the  electron  is  quickly  lost  again  to  dissociative 
attachment  to  F^.  A chain  length  of  3.5  is  proposed  in  Reference  20  for 
this  process.  This  reaction  was  not  included  in  the  simple  plasma  model 
of  the  previous  section*,  however,  if  it  were  included,  a larger  value  of 
the  chlorine  dissociative  attachment  rate  could  be  used  along  with  a corre- 
spondingly smaller  dissociation  fraction  by  the  discharge.  These  changes 
would  be  small,  less  than  factors  of  3 to  4 in  the  calculated  dissociative 
attachment  rate. 

B . Chemical  Kinetics  and  Laser  Gain  Computer  Model 

The  kinetic  model  included  only  the  first  four  HC1  vibrational  energy 
levels  (v=0,l,2,3).  Only  single  quantum  energy  exchange  and  deactivation 
processes  were  considered.  HCL  V-V  ladder  climbing  processes  were  limited 
to  HC 1 (v=l)  + HC1  (v')  X (v'+l)  + HC1  (0)  with  a partially  reflecting 
boundary  for  v'  = 3,  i.e.,  HC1  (v=4)  molecules  were  considered  lost  via  fur- 
ther V-V  processes.  The  vibrational ly  excited  formed  by  the  discharge 
was  assumed  to  be  in  a Boltzmann  equilibrium  and  only  (v=0,l)  were  con- 
sidered as  reacting  species. 

A recent  review  of  rate  constants  relevant  to  the  HC1  chemical  laser 

21 

system  by  N.  Cohen  and  J.  Bott  was  extremely  useful  in  summarizing  and 
providing  a critical  review  of  the  literature.  These  data  were  supplemented 
by  more  recent  data  appearing  after  the  publication  of  Reference  21.  Table 
3 summarizes  the  rate  coefficients  for  the  chemical  reactions  and  deactiva- 
tion processes  included  in  the  model. 


Rate  Coefficients  for  the  HC1  Chemical  Laser  System 
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The  kinetic  master  equation  included  all  the  reactant  and  product 
species  and  the  HC1  and  vibrational  level  populations  within  the  limita- 
tions discussed  above.  In  the  case  of  stimulated  emission,  there  exists 
an  additional  set  of  coupled  equations  for  the  HC 1 laser  transitions  which 
are  included  in  determining  the  vibrational  level  populations.  These  kinetic 
and  laser  equations  were  solved  simultaneously  with  the  energy  conservation 
equation  using  the  Runge-Kutta-Gi 1 1 method  of  integration. 


"*5  37 

Calculations  for  H“  Cl  and  H Cl  were  run  separately  using  a constant 

(0.75  or  0.25)  in  the  kinetics  model  which  represents  the  fractional  isotopic 
35  37 

abundance  of  Cl  or  Cl.  The  constant  multiplied  all  production  terms  for 
specific  vibrational  levels  in  HC1  in  such  a manner  that  the  total  HC1  con- 
centration still  appears  as  a deactivating  species  and  the  correct  chlorine 
atom  and  molecule  densities  were  computed. 

3b  37 

The  model  does  not  include  V-V  exchange  between  H Cl  and  H Cl.  The 

23 

exchange  rate  reported  by  Leone  and  Moore  (1.1  x 1 0 1 3 cm3/mole-sec)  implies  a 
V-V  time  of  10  to  100  psec  for  HC1  densities  on  the  order  of  10-9  to  10-8 
moles/cc.  Hence,  for  the  range  of  HC1  (v=0,l)  densities  found  in  the  model, 
the  rate  of  isotopic  V-V  exchange  is  slow  on  the  time  scale  of  interest. 

In  order  to  calculate  the  laser  flux  from  an  optical  cavity,  it  is 
necessary  to  know  the  gain  coefficient  of  the  amplifying  medium.  For  the 
experimental  condition  of  this  study,  the  gain  coefficient  is  fully  pressure 
broadened  and  is  given  for  the  P-branch  (J-l  -*•  J)  by 


UL  4tt 
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where 

= transition  line  center  wavelength 

A^  = spontaneous  transition  probability 

B.j ,B.  = characteristic  rotational  constants 
L for  upper  and  lower  levels 

nU,nL  = tota^  population  of  upper  and  lower 
vibrational  levels,  respectively 

Tr  = rotational  temperature 

Av  = Lorentz  broadened  linewidth  (full  width 
at  half  maximum) 

J = lower  level  rotational  quantum  number 

The  wavelengths  and  rotational  constants  were  obtained  from 

G.  Herzberg,  the  transition  probabilities  from  J.  M.  Herbelin  and  G. 

34  35 

Emanuel,  and  the  Lorentz  broadening  from  H.  Babrov  et  al . The  Lorentz 

broadening  of  HC1  is  strongly  J-dependent  for  not  only  self-broadening  but 

also  for  H 2 and  Ar  broadening.  Since  the  experimental  data  of  Reference 

36  for  and  Ar  did  not  cover  all  the  J values  of  interest,  the  broadening 

had  to  be  extrapolated  to  the  higher  J's  as  shown  in  Figure  14  . For  HC1 

self-broadening, the  experimental  data  of  Benedict  et  al.  was  used  for 

J < 13;  for  larger  values  of  J,  the  value  at  P ( 1 3)  was  used  since  this  was 

essentially  the  asymptote  of  Avc  as  a function  of  P(J). 

The  laser  flux  in  the  cavity  may  be  found  by  integrating  the  equation 
for  cavity  intensity  ($) 


= ct(g-gth)(LCAV/Lm) 

where  c is  the  speed  of  light  and  g**1  is  the  threshold  gain  (where  the 
cavity  gain  balances  the  losses  exactly).  The  final  term,  which  is  the  ratio 


Av  (HWHM)  (cnr Vatm) 
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.1  _ J P-BRANCH 
|m|  “ J+l  R- BRANCH 

Figure  14.  HCl-Hp  and  HCl-Ar  Optical  Broadening  Coefficients 
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of  the  active  medium  cavity  length  (50  cm)  to  the  mirror  separation  length 
(130  cm),  accounts  for  the  fraction  of  time  the  laser  mode  spends  in  the 
active  region.  In  the  case  of  HC1,  where  several  vibrational  lines  lase 
simultaneously,  a cavity  intensity  equation  was  used  for  each  vibrational 
line.  In  addition,  each  vibrational  line  would  have  an  equation  for  each 
rotational  line;  however,  this  was  simplified  by  solving  only  one  cavity 
intensity  equation  for  each  vibrational  line,  with  the  transition  going  on 
the  optimum  gain,  rotational-vibrational  transition  (i.e.,  complete  rota- 
tional equilibrium). 

C.  Computer  Model  Comparisons  with  Experimental  Observations 

The  electron  beam-initiated  chemical  laser  presents  the  least  number 
of  variables  for  modelling  the  chemical  kinetics.  The  principal  V-T  rates 
for  Cl  and  H atoms  on  HCl(v)  have  been  measured  by  several  independent 
groups,  and  the  rates  are  in  good  agreement.  The  V-V  rates,  at  least  for 
v = 1,  for  HCl(v)  and  (^(v)  are  also  well  established  and  were  not  altered 
in  the  present  calculations. 

Four  rates  were  found  to  be  extremely  sensitive  in  analyzing  the  elec- 
tron beam-initiated  chemical  laser.  Production  rates  of  H and  Cl  atoms  and 
the  rate  constants  for  the  "hot"  H+Clj,  and  "cold"  Cl+H^  reactions  were 
investigated  in  some  detail. 

The  maximum  amount  of  H-atom  production  was  bounded  by  the  energy 
available  from  the  primary  electron  beam  (~30  eV/ion  pair),  which  implies  a 
maximum  of  six  hydrogen  atoms  (4.5  eV/H  atom)  per  ion  pair.  The  lower 
bound  on  the  H-atom  production  is  established  by  the  observed  laser  energy 
which  was  almost  directly  proportional  to  the  H-atom  production  rate  and  the 
rate  constant  chosen  for  the  H+C^  reaction.  However,  as  the  reaction  rate 
was  increased  in  the  model  for  decreased  H-atom  production,  the  overall 
pulse  shape  changed.  In  addition,  it  was  found  that  good  agreement  of  the 
model  with  the  discharge  experiments  could  only  be  obtained  by  maintaining 
the  H+Cl^  rate  at  the  value  shown  in  Table  3 . This  rate  coefficient  is 
a factor  of  1.75  times  larger  than  the  value  recommended  by  Reference  21 
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and  well  within  the  overall  spread  of  experimental  results  summarized  in 
Reference  21.  The  H-atom  production  rate,  which  was  consistent  with  this 
choice  of  rate  constant,  was  2.25  H-atoms  per  ion  pair.  Reference  20  esti- 
mates (Eq.  A19)  that  two  H-atoms  would  be  produced  per  ion  pair,  which  is 
in  good  agreement  with  the  value  derived  in  the  present  work. 
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Laser  emission  from  2-1  H Cl  also  established  a lower  bound  to  the 

amount  of  H-atom  present,  since  its  time  to  threshold  was  quite  sensitive 

to  the  H-atom  production.  If  the  H-atom  production  was  decreased  to  one 

37 

H-atom  per  ion  pair,  the  calculated  2-1  H Cl  threshold  time  was  delayed 
too  long  (0.8  psec)  from  the  experimentally  observed  time. 

The  cold  reaction  (Cl+H^)  rate  constant  listed  in  Table  3 is  a fac- 
tor of  2 larger  than  the  recommended  value  of  Reference  21.  However,  a 
38 

recent  reevaluation  of  this  rate  constant  is  within  20  percent  of  the  value 
used  in  the  present  work.  The  magnitude  of  this  rate  coefficient  controls 
the  1-0  laser  threshold  time  since  the  reaction  produces  only  HC1(0).  If 
the  production  rate  of  HC 1(0)  by  the  cold  reaction  was  reduced,  the  1-0 
laser  transitions  reached  threshold  too  early  and  the  pattern  of  three 
prominent  peaks  (Fig.  4(a))  was  changed  to  only  two. 

For  the  electron  beam-initiated  chemical  laser,  the  Cl  atom  produc- 
tion rate  by  dissociative  attachment  would  be  limited  to  the  ion-pair 
production  rate,  whereas  it  was  found  that  a rate  four  times  higher  was 
required  to  be  consistent  with  both  the  experimentally  determined  burn  rate 
of  C 1 p and  the  kinetics  of  the  HC1  laser.  Therefore,  it  was  assumed  that 
processes  similar  to  the  H^  dissociation  by  the  e-beam  were  responsible  for 
producing  an  additional  four  chlorine  atoms  per  ion  pair.  Comparisons  of 
the  measured  C 1 2 disappearance  rate  and  the  calculated  Cl2  time-dependent 
concentration  are  shown  in  Figure  15.  The  measured  rate  of  C 1 ^ consump- 
tion, 1 to  2 x 10“ 3 moles/c.c.-sec,  was  sensitive  to  the  electron  beam 
pulse  length  as  shown  in  Figure  15.  Independent  measurements  of  the  primary 
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Figure  15.  Rate  of  Chlorine  Disappearance  Using  Only  the  Electron  Beam 
for  Excitation.  Gas  pressure  was  400  torr  and  the  gas  mix- 
ture was  Ar/f^/C^  (49/49/2). 
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electron  beam  current  density  time  behavior  showed  that  the  profile  is  some- 
what skewed  so  that  the  maximum  current  density  occurs  late  in  the  pulse. 

The  computer  model  assumes  a square  pulse  current  density  time  behavior 
and  will  therefore  predict  only  small  changes  in  the  burn  rate  for 
correspondingly  small  changes  in  the  electron  beam  pulse  duration.  The 
agreement  between  the  measured  and  calculated  consumption  rate  was  rea- 
sonable based  on  the  above  assumption  and  the  value  chosen  for  the  cold 
reaction  rate  constant. 

The  additional  chlorine  atom  production  is  nearly  identical  to  that 

20  39 

found  in  two  independent  studies  ’ of  F-atom  production  by  direct  e-beam 

irradiation  of  dilute  F„  rare  gas  mixtures.  Penning  ionization  and  charge 

^ + 

transfer  mechanisms  were  proposed  in  Reference  40  for  F?  production  by 
+ + ^ 

He*,  He~  , and  Ar  , followed  by  ion-ion  and  ion-electron  recombination  re- 
^ + 

actions  of  F^  which  lead  to  the  production  of  F-atoms. 

A calculated  laser  pulse  time  history  is  shown  in  Figure  16  for 
electron  beam  initiation  of  the  HC1  chemical  laser.  The  agreement  with  ex- 
periment is  quite  good  (see  Fig.  4(a))  and  the  calculated  laser  pulse  energy 

(20  mJ)  is  well  within  the  range  of  experimentally  measured  energies  (Table  2). 
37 

The  calculated  H Cl  2-1  laser  threshold  time  and  intensity  are  also  shown 
and  are  in  excellent  agreement  with  the  data  presented  in  Figure  5. 

35 

There  typically  was  a weak,  late  H Cl  1-0  P(9)  feature  in  the 
electron  beam-initiated  laser  pulse  occurring  1 psec  after  the  last  strong 
spike  (Fig.  4(a)).  Efforts  to  obtain  this  feature  in  the  computer-simulated 
pulse  were  unsuccessful.  A wide  range  of  H-atom  production  rates  and 
H+Cl ^ reaction  rate  constants  were  tried.  Nonuniformity  of  the  electron 
beam  spatial  distribution  was  a possible  explanation  for  the  late  line. 
Calculations  simulating  a region  of  lower  e-beam  current  density  showed 
that  essentially  only  1-0  transitions  lase  at  approximately  2 psec  after 
the  termination  of  the  electron  beam  pulse,  which  is  the  time  at  which  the 
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experimentally  observed  1-0  P (9 ) line  appears.  However,  experimental  mea- 
surements of  the  electron  beam  density  profile  showed  that  spatial  nonuni- 
formities would  not  have  been  significant  across  the  area  of  interest. 

Hence,  the  observed  late  1-0  P(9)  is  not  suitably  explained  by  nonuniformities 

0 'll 

transverse  to  the  optical  axis.  However,  as  indicated  by  Searles,  non- 
uniformities  can  occur  as  gradients  associated  with  the  electron  beam  energy 
deposition  in  the  laser  gas.  This  type  of  nonuniformity  is  not  readily 
measured  and  must  be  evaluated  numerically. 

Values  for  the  vibrationally  enhanced  reaction  (Cl+H2(v))  rate  coeffi- 
cient and  the  production  rate  of  H and  Cl  atoms  by  the  discharge  were  re- 
quired for  the  analysis  of  the  electron  beam-stabilized  discharge  experiments. 

Initial  values  for  the  enhanced  reaction  rate  constant  were  taken  from 
5 

Wilkins'  recent  calculations.  The  H-atom  and  H2(v)  production  rates  were 
determined  by  solutions  to  the  electron  Boltzmann  equation  at  various  values 
of  E/N.  The  Cl -a tom  production  rate  by  the  discharge  was  found  by  parametric 
considerations  of  its  effect  on  the  laser  pulse  shape,  laser  energy,  and 
disappearance  rate  of  the  Cl.,. 

It  was  immediately  obvious,  as  shown  in  Figure  17,  that  the  Cl+H2(v) 
rate  constant  derived  in  Reference  5 was  much  too  large.  The  calculated 
laser  pulse  lasts  6.5  psec  whereas,  in  fact,  it  was  3.5  psec  long.  The  cal- 
culated output  energy  (600  mJ)  was  also  much  too  large.  Including  chlorine 
dissociation  in  the  model  (Fig.  18)  reduced  the  pulse  duration  since  the  Cl 
atoms  have  rapid  V-T  deactivation  rates  on  HCl(v).  However,  the  isotopic 
H Cl  laser  threshold  times  were  much  too  early  (0.7  psec)  compared  to  the 
experimentally  measured  values  presented  in  Figure  6,  and  the  laser  energy 
was  again  much  too  large  (750  mJ).  A limit  on  the  extent  of  chlorine  dissoci- 
ation was  the  measured  disappearance  rate  of  Cl2  and  the  initial  absolute 
decrease  in  the  Cl2  concentration  immediately  following  the  discharge.  This 
limit  on  the  chlorine-atom  production  rate  forced  the  Cl+H2(v)  rate  to  be 
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TIME  (usee) 

Figure  17.  Calculated  HC1  Laser  Pulse  Shape  for  Low  E/N  (1.55  x 10“16 
V- cin2)  Discharge  Excitation.  The  rate  constant  for  the 
vibrational ly  enhanced  reaction  was  taken  from  Reference 
5,  and  no  chlorine  atom  production  by  the  discharge  was 
assumed. 
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decreased  by  a factor  of  10  below  that  calculated  in  Reference  5,  and  the 
resultant  laser  pulse  shape  is  shown  in  Figure  18.  If  the  rate  constant 
for  the  reaction  was  taken  to  be  100  times  smaller  than  the  value  in  Reference 
5 and  chlorine  dissociation  limited  to  dissociative  attachment,  then  the  laser 
pulse  merely  consists  of  a single  sharp  spike  with  the  total  energy  wholly 
determined  by  the  H-atom  production  rate. 

Figures  19  and  20  show  the  calculated  laser  pulse  shapes  using  1/10 
the  calculated  value  of  Reference  5 for  k^(v)  'for  the  low  and  high  E/N  cases 
presented  in  Figure  4(b)  and  (c).  The  overall  agreement  is  quite  good,  and 
the  calculated  laser  energies  are  close  to  the  range  of  experimentally  mea- 
sured values  shown  in  Table -2.  In  Figure  21  the  experimental  time-resolved 

35 

laser  intensity  profiles  for  the  H Cl  1-0,  2-1,  and  3-2  bands  are  plotted 
along  with  the  calculated  profiles  for  the  low  E/N  discharge  case.  As  seen 

OC 

in  Figure  21(b),  the  calculated  1-0  H Cl  laser  transition  intensity  is 

35 

approximately  twice  that  for  the  2-1  H Cl  band.  Whereas  experimentally, 
it  was  found  that  the  2-1  transition  was  slightly  more  intense  early  in  the 
laser  pulse.  A similar  contradiction  is  found  for  the  high  E/N  discharge 
case  as  well.  Additional  calculations  have  shown  that  the  assumption  of  a 
large  HCl(v=l)  product  fraction  in  the  Cl+I^v)  reaction  is  not  the  cause  of 
the  disagreement.  The  discrepancy  could  not  be  remedied  by  assuming  only 
an  8 percent  fractional  energy  distribution  into  v = 1 from  the  hot  reaction 
(see  Table  3).  The  experimental  summation  shown  in  Figure  21  is  quite 
sensitive  to  the  exact  peak  intensities  and  time  histories  for  all  the  lines. 
Hence,  a slight  shift  (100  to  200  nsec)  in  a line  position  could  strongly 
affect  the  total  band  intensity  profile,  yet  only  slightly  alter  a complete 
composite  reconstruction  as  presented  in  Figure  8.  Therefore,  at  present, 
the  discrepancy  presented  in  Figure  21  is  assigned  to  an  artifact  of  the 
data  reconstruction  procedure. 

The  calculated  time  history  of  a particular  P(J)  transition  was 
typically  in  error  by  + 1 or  +2  in  J at  a given  time.  This  was  not  too 
surprising  considering  the  sensitivity  of  the  gain  on  a particular  J to 
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Figure  18.  Calculated  HC1  Laser  Pulse  Shape  for  Low  E/N  (1.55  x 10~16  V-cm  ) 
Discharge  Excitation.  Chlorine  atom  production  by  the  discharge 
was  10  percent  of  the  discharge  energy  input  for  both  cases  shown. 


77  00752 


LASER  INTENSITY  (<5>  x 105  watts/cm 


44 


Figure  19.  Calculated  HC1  Chemical  Laser  Pulse  Shape  for  Low  E/N 
(1.55  x 10-16  V-cm2)  Discharge  Excitation  (Current 
Density,  24  A/cm2) 
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Figure  20.  Calculated  HC1  Chemical  Laser  Pulse  Shape  for  High  E/N 
(2.5  x 10-16  V-cm2)  Discharge  Excitation  (Current 
Density,  24  A/cm2) 
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Figure  21.  Comparison  of  Experimentally  Determined  Band  Intensities  with 
Calculated  Profiles  for  the  1-0,  2-1,  and  3-2  H^Sci  Laser 
Transitions  for  the  Low  E/N  (1.55  x 10-16  V-cm2)  Discharge  Case 
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the  inversion  ratio  of  the  two  vibrational  levels  involved.  In  addition, 
it  was  found  that  the  J time  sequence  could  be  strongly  altered,  depending 
on  the  choice  of  the  collision  broadening  parameter.  Further  refinements 
in  the  model  would  be  necessary  for  a more  accurate  tracking  of  the  de- 
tailed P(J)  time  history. 

37 

The  laser  threshold  times  of  the  H Cl  isotopes  provided  very  sensi- 
tive criteria  for  establishing  the  H-atom  production  rates  by  the  discharge. 
For  example,  in  the  low  E/N  calculations  (Fig.  19),  the  2-1  H Cl  laser 
transitions  did  not  reach  threshold  unless  the  H-atom  production  exceeded 

5 percent  of  the  discharge  energy,  in  addition  to  the  H-atom  production  by 
the  electron  beam  itself.  This  production  rate  by  the  discharge  essentially 
equals  the  electron  beam  H-atom  production.  It  could  be  argued  that  the  de- 
activation rates  chosen  for  HCl(v=2)  were  too  fast  and  therefore  higher 
H-atom  production  was  necessary  to  overcome  the  losses.  However,  if  the 

deactivation  rates  for  v = 2 were  decreased,  via  lower  chlorine  atom  produc- 
35 

tion,the  H Cl  2-1  laser  transitions  did  not  turn  off  by  2.2  usee  as  found 

35 

in  the  experimental  time  history  (Fig.  5).  The  calculated  H Cl  2-1  laser 
intensity  drops  to  1/10  its  maximum  by  3.0  ysec,  which  indicates  that  the 
deactivation  processes  for  v = 2 were,  in  fact,  slightly  underestimated  and 
therefore  the  H-atom  production  rate  may  also  have  been  underestimated. 

Hydrogen  dissociation  by  the  low  E/N  discharge  was  calculated  (see 
Sec.  IV-A)  to  be  0.1  percent  of  the  total  discharge  energy  input,  although 

6 percent  of  the  discharge  energy  needed  to  be  used  for  hydrogen  dissocia- 
tion as  discussed  above.  The  high  E/N  Boltzmann  calculations  agreed  very 
well  with  the  9 percent  discharge  energy  fraction  into  dissociation  used 

in  the  kinetics  model.  The  large  discrepancy  at  low  E/N  cannot  be  explained 
at  present. 
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The  chlorine  disappearance  rate  helped  to  constrain  the  vibrational ly 
enhanced  rate  constant.  A simple  kinetic  analysis  using  the  first  three 
reactions  in  Table  1 yields  for  the  C 1 2 burn  rate: 

dCl?  r -| 

~dT  = [kH?(v)^H2^v^+kcold^H2 J 

assuming  that  the  chlorine  atom  density  reaches  a steady  state  value  follow- 
ing the  discharge  pulse.  An  estimate  for  the  chlorine  atom  steady  state 
concentration,  as  verified  in  the  computer  model  calculations,  can  be  ob- 
tained from  Figure  22  using  the  initial  drop  in  C 1 2 density  at  approximately 
1 ysec,  which  leads  to  a value  of  4 x 10"8  moles/c.c.  The  vibrational 
temperature  of  Hg  calculated  immediately  after  the  discharge  is  approximately 
2000  °K,  hence  the  H2(v)/H2(0)  ratio  is  0.05.  Therefore,  the  experimentally 
measured  Cl2  disappearance  rate  (E/N  = 1 . 6 x 1 0“ 1 6 V-cm2)  of  3.4  x 10~2 
moles/c.c. -sec,  implies  a rate  ratio  (^(v^H  (0)  12-5-  Note  that  this 

is  the  rate  enhancement  at  350  °K,  which  is  the  calculated  temperature  at 
the  end  of  the  discharge  pulse.  A rate  ratio  of  20  at  350  °K  can  be 
calculated  from  the  rate  constants  in  Table  3.  This  is  reasonably  good 
agreement  between  the  rate  evaluation  based  on  laser  characteristics  and 
that  based  on  a direct  C 1 2 removal  rate  measurement. 

The  calculated  chlorine  disappearance  rate  using  the  complete 
kinetics  code  was  4 x 10'2  moles/c.c. -sec,  which  is  in  very  good  agreement 
with  the  measured  values,  which  ranged  from  3.4  x 10'2  to  3.9  x 10-2  moles/ 
c.c.-sec.  For  the  high  E/N  case,  the  measured  disappearance  rate  was 
approximately  9 x 10"2  moles/c.c. -sec  and  the  calculated  rate  was  11.5  x 10~2 
moles/c.c. -sec.  Ten  percent  of  the  discharge  energy  input  for  the  high  and 
low  E/N  cases  was  placed  in  chlorine  dissociation.  For  the  low  E/N  case, 
this  corresponds  to  4 percent  direct  chlorine  dissociation  by  the  discharge. 

5 

Wilkins'  prediction  of  a vibrationally  excited  product  molecule  in 
the  reaction  Cl+H2(v)  -*•  MCI  (0,1)  + H appears  to  be  correct.  If  the  rate  of 
the  vibrationally  enhanced  reaction  does  not  include  HC 1(1)  production,  the 
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late  time  1-0  laser  transitions  (see  Fig.  6)  do  not  occur.  Wilkins'  cal- 
culations predict  approximately  a ratio  of  3 for  the  rate  constant  producing 
HC 1(1)  to  that  for  HC1(0),  which  seems  consistent  with  the  experimental 
observations  and  computer  model. 

A sink  term  for  the  HCl(v)  levels  was  introduced  above  v = 3.  Since 
the  chemically-produced  contribution  to  HCl(v)  for  v > 3 was  less  than  5 
percent,  and  these  levels  do  not  lase;  these  levels  were  not  included  in  the 
calculations  as  discussed  earlier.  A sink  was  introduced  into  the  calcula- 
tions for  the  V-V  ladder  climbing  reaction:  HC 1(1)  + HCL(3)  J HC1 (4)  + HC1(0), 

and  HC 1 (4)  molecules  were  then  assumed  not  to  reenter  the  laser  kinetics. 

O 

The  rate  constant  shown  in  Table  3 represents  one-third  of  the  Vc  scaling 
law  commonly  used  and  gave  the  best  agreement  with  the  observed  duration 
of  the  1-0  laser  transitions  at  the  end  of  the  laser  pulse. 


SECTION  V 


CONCLUSIONS 

Experimental  measurements  of  an  electrically-initiated  HC1  chemical 

laser  system  have  been  compared  with  a reaction  kinetics-laser  gain  computer 

model.  An  estimate  for  the  vi brational ly-enhanced  rate  coefficient  for  the 

reaction  Cl+H9(l)  -*■  HC1(0,1)  + H has  been  derived.  It  is  ten  times  smaller 
£ 5 

than  a recent  theoretical  prediction.  However,  the  rate  constant  found 

is  40  times  larger  than  the  C 1 +H2 (0)  rate  given  in  Table  1 (T=300  °K). 

This  finding  is  in  qualitative  agreement  with  an  earlier  experimental  mea- 

41 

surement  of  Stedman  et  al.  • where  a vibrational  rate  enhancement  of  "two 

orders  of  magnitude"  was  reported  in  a mass  spectrometric  discharge  flow 

3 

experiment.  Sims  et  al.  report  a theoretical  rate  enhancement  of  23  times 

(T=300  °K),  which  is  only  a factor  of  2 smaller  than  the  present  findings. 

4 

Persky  and  Baer  calculated  a rate  enhancement  of  35  times,  encouragingly 
close  to  the  present  evaluation. 

Since  the  C 1 +H2 (0 ) rate  constant  is  strongly  temperature  dependent 
(see  Table  3),  the  relative  rate  enhancement  k(H2(v))/k(H2(0))  decreases 
rapidly  with  increasing  temperature,  for  example  at  400  °K,  the  rate  ratio 

will  decrease  to  9,  using  Wilkins'^  temperature  dependence  for  the  Cl+H9(v) 

4 c 

reaction.  Persky  and  Baer  calculate  a temperature-dependent  rate  constant 

3 

ratio  at  400  °K  of  13,  whilst  Sims  et  al.  report  a ratio  of  12. 

Modelling  of  experiments  using  electron  beam  only  initiation  of  the 
chemical  laser  system  was  dependent  upon  the  Cl  and  H-atom  production 
rates  and  the  H+Cl^  and  Cl+H^  reaction  rates.  The  reaction  rate  coefficients, 
which  gave  the  best  agreement  with  experimental  measurements,  are  shown  in 
Table  3.  Estimates  for  the  chlorine  dissociation  fraction  in  both  the 
electron  beam  and  discharge-excited  gas  mixtures  were  based  on  measurements 
of  the  chlorine  disappearance  rate.  Hydrogen-atom  production  by  the  elec- 
tron beam  was  directly  related  to  the  laser  energy  and  pulse  shape  and  could 
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be  fitted  to  the  experimental  observations  only  over  a fairly  narrow 

37 

range  (less  than  a factor  of  2).  The  H Cl  laser  threshold  time  also  pro- 
vided a valuable  benchmark  for  the  H-atom  production  rate.  Hydrogen-atom 
production  in  the  discharge  experiments  were  estimated  a priori  using  pub- 
lished cross  sections  for  hydrogen  excitation  channels  in  a Boltzmann  cal- 
culation. Good  agreement  with  the  observed  HC1  laser  energy  and  spectro- 
scopy was  obtained  only  by  increasing  the  calculated  dissociation  fraction 
by  a factor  of  60  for  the  low  E/N  cases.  For  the  high  E/N  discharge  ex- 
periments, the  proposed  dissociation  fraction  was  essentially  the  same  as 
the  value  found  in  the  Boltzmann  calculation. 

The  value  obtained  in  this  study  for  the  vibrationally-enhanced  rate 
coefficient  was  constrained  by  the  above  considerations  for  the  H+Cl^  and 
Cl+H^  rate  coefficients  and  the  H and  Cl-atom  production  rates.  It  was 
found  that  the  rate  constant  fitted  to  the  experimental  results  could  not 
be  varied  by  more  than  20  to  30  percent  without  serious  discrepancies  appear- 

5 

ing  in  the  laser  pulse  shape  comparison.  Wilkins’  theoretical  prediction 
for  the  formation  of  HC 1(1)  in  the  Cl+H£(v)  reaction  is  consistent  with 
experimental  observation  of  a long"tail"  on  the  laser  pulse,  consisting 
mainly  of  1-0  laser  transitions  which  could  be  duplicated  in  the  present 
calculations  of  the  laser  pulse  time  history. 

An  immediate  practical  consequence  of  the  smaller  rate  constant 
k^  ^ is  the  inherent  inefficiency  of  the  electrically-excited  HC1  chemical 
laler.  Most  of  the  energy  stored  in  the  H?(v)  is  never  recovered  on  the 
time  scale  of  the  laser  pulse.  The  deactivation  rate  of  Cl  atoms  predomi- 
nates over  other  deactivation  routes  under  the  conditions  of  these  experi- 
ments and  determines  the  overall  time  constant  for  the  laser  pulse  length. 

The  maximum  electrical  to  laser  energy  conversion  efficiency  obtained  was 
4.5  percent. 
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SECTION  I 


INTRODUCTION 

Self-sustained  and  UV  preionized  discharge  techniques  for  excita- 
tion of  the  rare  gas  halide  lasers,  e.g.,  KrF,  XeF,  and  ArF,  have  proved 
very  useful  for  small  laser  systems  but  appear  to  be  limited  to  pulse 
energies  less  than  1 J.  At  present,  the  only  viable  scheme  for  higher 
pulse  energies  and  high  average  power  appears  to  be  the  electron-beam  sus- 
tained discharge.  It  has  been  shown  that  the  discharge  physics  in  mixtures 
of  the  rare  gases  and  halogen  containing  molecules  is  dominated  by  electron 
impact  excitation  and  ionization  of  the  rare  gas  metastables,  leading  to  an 
effective  upper  bound  on  the  applied  E/N  of  the  discharge  to  avoid  instabil- 
ities.^ In  order  to  achieve  long-term  stable  discharges,  it  appears  necessary 
to  operate  at  modest  pump  power  density  and  a sufficiently  large  gain  length 
in  order  to  saturate  the  medium.  This  is  in  contrast  with  the  early  experi- 
ments on  electron-beam  ionized  discharge  excitation  of  KrF  which  involved 
input  power  densities  of  the  order  of  1 MW/cm3  and  excitation  pulse  length 
limited  to  approximately  100  nsec  by  discharge  instabilities.  The  present 
measurements  were  undertaken  in  order  to  investigate  the  minimum  threshold 
power  density  required  to  achieve  stimulated  emission  in  both  KrF  and  XeF 
for  long  pulse  excitation  (~  1 psec)  by  the  electron-beam  sustained  dis- 
charge technique. 


SECTION  II 


DESCRIPTION  OF  EXPERIMENT 


A schematic  of  the  experimental  apparatus  is  shown  in  Figure  1.  The 
cold  cathode  electron  gun  was  described  in  the  previous  section  dealing  with 
the  HC1  laser  work.  It  was  operated  with  a 3 cm  cathode  to  anode  spacing  at 
a Marx  voltage  of  160  kV.  This  resulted  in  a 1-psec,  4-amp/cm2  ionization 
pulse  for  the  discharge.  The  current  density  through  the  aluminized  Kapton 
foil  was  estimated  by  simul taneously  measuring  the  energy  deposition  with 
a calorimeter  and  the  temporal  history  with  a current  transformer.  The 
e-beam  uniformity  was  verified  by  firing  the  gun  into  a piece  of  Ozalid 
paper  placed  adjacent  to  the  discharge  cathode  screen.  It  was  found  nec- 
essary to  sandpaper  the  graphite  cathode  on  the  electron  gun  each  time  a 
foil  ruptured  in  order  to  ensure  a uniform  current  distribution  through 
the  foil . 

The  discharge  chamber  was  constructed  of  high  density  polyethylene, 
which  provides  a passive  surface  for  the  highly  reactive  fluorine.  The 
nominal  discharge  volume  is  625  cm3  and  is  formed  by  a 50-cm  x 5-cm  alumi- 
num anode  spaced  2.5  cm  from  a screen  cathode  adjacent  to  the  electron  beam 
foil  window.  The  discharge  cell  was  fitted  with  fused  silica  or  calcium 
fluoride  windows  at  Brewster's  angle,  leaving  an  inactive  12  cm  region  at 
each  end  of  the  discharge.  The  fused  silica  windows  were  attacked  by  the 
fluorine  containing  gas  mixtures  and  had  to  be  changed  frequently,  while 
the  calcium  fluoride  windows  showed  little  evidence  of  damage.  A stable 
optical  cavity  was  obtained  by  the  use  of  2-m  radius-of-curvature  mirrors 
spaced  1.3  m apart  or  a 2-m  radius-of-curvature  mirror  opposite  a flat 
mirror  at  the  same  separation.  The  maximum  reflectance  mirrors  for  the 

O 

KrF  experiments  (coated  for  2500  A)  had  a loss  of  3 percent  while  the  cor- 

O 

responding  mirrors  for  XeF  (coated  for  3500  A)  had  less  than  1 percent  loss. 
Output  coupling  mirrors  with  transmissions  of  12  and  30  percent  were  also 
available  for  the  KrF  and  XeF  experiments,  respectively. 
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Figure  1.  Experimental  Layout  for  KrF  and  XeF  Threshold  Pumping  Measurements 
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All  of  the  present  experiments  were  performed  at  a total  gas  pressure 
of  1 atm  in  order  to  minimize  the  number  of  foil  failures.  The  KrF  experi- 
ments were  performed  in  mixtures  of  Kr  and  diluted  in  Ar  while  most  of  the 
XeF  experiments  were  performed  in  mixtures  of  Xe  and  NF^  diluted  in  Ar.  Some 
preliminary  experiments  using  NF3  and  F^  as  the  F donor  for  the  XeF  experi- 
ments indicated  higher  efficiency  with  NF3.  The  Ar  diluent  was  used  to  pro- 
vide a high  stopping  power  for  the  electron  beam  ionization. 

The  gases  were  prepared  in  a high-purity  stainless  steel  manifold 
and  then  allowed  to  sit  until  diffusive  mixing  had  taken  place.  The  en- 
tire manifold  was  thoroughly  passivated  by  filling  it  with  gradually  in- 
creasing fluorine  concentrations  ending  up  with  an  exposure  of  several 
hours  at  a fluorine  pressure  of  250  torr.  Disposal  of  the  fluorine  was 
accomplished  by  pumping  it  slowly  through  a charcoal  trap  where  it  rapidly 
reacts  with  the  charcoal  to  form  volatile  fluorinated  hydrocarbons.  The 
discharge  chamber  was  passivated  by  following  the  same  procedure  described 
for  the  gas  handling  manifold.  The  effectiveness  of  this  passivation  was 
verified  by  placing  an  Ar  + 1%  F,  gas  mixture  in  the  discharge  chamber  and 
monitoring  the  absorption  of  the  Hg  2537  A resonance  line  by  F^.  Before 
the  discharge  chamber  was  passivated,  the  F^  disappeared  within  two  or 
three  minutes  after  being  added.  However,  after  the  discharge  chamber  was 
cleaned  and  passivated,  there  was  no  detectable  loss  of  F^  for  times  as 
long  as  20  minutes  after  filling,  even  when  several  e-beam  sustained  discharge 
pulses  were  fired  into  the  gas.  The  UV  absorption  scheme  was  also  used  to 
verify  the  composition  of  the  premixed  gases  and  indicated  that  the  gas  mix- 
ing manifold  was  well  passivated. 

The  discharge  energy  for  these  experiments  was  supplied  with  two 
different  pulse  forming  networks--one  with  a pulse  duration  of  approximately 
100  nsec  and  the  other  with  a pulse  duration  of  about  1 psec.  The  initial 
work  was  carried  out  with  the  100  nsec  network  because  it  was  possible  to 
construct  the  pulser  from  parts  already  on  hand  while  we  waited  for  delivery 
of  the  capacitors  for  the  1 psec  line.  The  100  nsec  network  was  constructed 
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from  twenty  2.4  nF  door-knob  capacitors  with  a voltage  rating  of  40  kV.  The 
capacitors  were  connected  between  two  parallel  plates  of  0.030  in.  thick  gal- 
vanized steel  and  were  arranged  in  five  rows  such  that  the  characteristic 
impedance  Z = /L/C-  ~ 1 SI  . The  1 psec  pulser  consisted  of  five  0.2  pf  capac- 
itors equally  spaced  along  a five-turn  solenoid  with  a coaxial  return  path, 
which  resulted  in  a characteristic  impedance  of  approximately  1 Si.  For  the 
threshold  measurements , two  of  the  capacitors  were  disconnected  from  one  end 
of  the  line,  which  shortened  the  pulse  duration  to  about  600  nsec.  The 
discharge  current  was  monitored  by  measuring  the  voltage  drop  across  the  low 
inductance  stainless  steel  resistor  in  series  with  the  return  side  of  the 
line.  The  discharge  voltage  was  measured  with  a resistive  divider  connected 
to  the  anode.  Comparison  of  this  divider  with  a Tektronix  (P6015)  high- 
voltage  probe  indicates-  that  it  had  a rise  time  of  50  nsec  or  less,  which  is 
adequate  for  the  pulse  durations  of  this  work.  The  pulse-forming  network 
was  discharged  into  the  gas  via  a spark  gap  triggered  at  a controlled  delay 
time  (usually  about  200  nsec)  following  the  initiation  of  the  e-beam  pulse. 


Voltage  and  current  monitors  were  used  on  the  electron  beam  and  the 
electric  discharge  while  a photodiode  (ITT  F4000  S-5)  was  used  as  the  primary 
diagnostic  to  monitor  the  cavity  flux  scattered  from  one  of  the  Brewster 
angle  windows.  Oscillograms  displaying  the  typical  temporal  variation  of 
these  parameters  are  shown  in  Figure  2.  The  threshold  pump  power  density 
was  determined  by  reducing  the  discharge  electric  field  to  the  point  where 
stimulated  emission  was  barely  discernible  within  the  discharge  pulse  width. 
One  drawback  in  these  experiments  was  the  necessity  to  vary  the  applied  E/N 
in  order  to  vary  the  input  power  density  at  the  fixed  current  density  of  the 
electron  gun.  The  applied  electric  fields  were  typically  in  the  range  of 
1.5  to  3 kV/cm-atm,  with  the  higher  electric  fields  being  achieved  at  the 
higher  concentrations  of  F£  and  NF^. 
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Figure  2.  Temporal  Variation  of  e-Beam  Voltage  and  Current,  Discharge 
Voltage  and  Current,  and  Photodiode  Signal  for  An  Ar  + 0.5% 
Xe  + 6.1%  NF3  Gas  Mixture  at  1 atm 


7 6 00302 


SECTION  III 
EXPERIMENTAL  RESULTS 


A . XeF  Laser  Threshold  Measurements 

Threshold  discharge  power  density  measurements  for  lasing  in  XeF  are 
shown  in  Figures  3a  and  3b  for  varying  NF^  concentration  at  fixed  Xe/NF^ 
ratio  and  for  varying  Xe/NF^  ratio  at  fixed  NF3  concentration.  Measurements 
were  made  with  both  the  maximum  reflectance  cavity  and  with  30  percent  out- 
put coupling.  Threshold  lasing  by  e-beam  excitation  alone  was  achieved  in 
the  maximum  reflectance  cavity  for  the  0.05%  and  0.1%  NF^  concentrations, 
the  e-beam  input  power  density  being  approximately  30  kW/cm3.  The  dis- 
charge threshold  measurements  for  these  cases  give  an  upper  bound  and 
correspond  to  the  power  density  necessary  to  significantly  increase  the 
output  and  to  obtain  the  threshold  delay  time  of  several  hundred  nanoseconds 
characteristic  of  all  the  measurements.  The  lack  of  lasing  with  e-beam  only 
excitation  for  NF,  concentrations  greater  than  0.1  percent  suggests  collis- 
ional  quenching  of  XeF  or  some  other  competition  involving  NF^  that  re- 
duces the  gain.  However,  in  no  case  was  threshold  lasing  achieved  by  the 
e-beam  alone  when  a 30  percent  output  coupler  was  used.  These  measurements, 
therefore,  represent  true  threshold  conditions  for  the  given  electron-beam 
current  density  and  excitation  pulse  length  of  the  order  of  1-psec.  The 
discharge  threshold  power  densities  are  typically  one-half  to  two-thirds  the 
e-beam  input  power  density. 

B . KrF  Laser  Threshold  Measurements 

Threshold  data  for  experiments  in  KrF  are  shown  in  Figure  4 for  F 2 
concentrations  ranging  from  0.1  to  0.8  mole  percent,  with  a fixed  Kr  con- 
centration of  6 percent.  These  mixtures  span  typical  gas  mixtures  used  in 
previous  e-beam  sustained  discharge  laser  experiments.  Measurements  were 
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Figure  3a.  Threshold  Electric  Discharge  Power  Density  for  XeF  Laser 
Emission  as  a Function  of  NF,  Concentration  for  Fixed 
Xe/NF3 
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Figure  3b.  Threshold  Electric  Discharge  Power  Density  for  XeF  Lasing 
as  a Function  of  Xe  Concentration  for  Fixed  NF^ 
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Figure  4.  Threshold  Electric  Discharge  Power  Density  Measurements 
for  KrF  Laser  Emission  as  a Function  of  F~  Concentration. 
The  Kr  concentration  was  fixed  at  6 percent. 
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made  with  both  the  maximum  reflectance  cavity  and  with  a 12  percent  output 
coupler.  Threshold  lasing  was  readily  achieved  by  the  e-beam  alone  for  F^ 
concentrations  of  0.2  percent  or  less  in  both  cavities  and,  therefore,  only 
an  upper  bound  can  be  set  to  the  threshold  discharge  power  density.  Dis- 
charge pumping  was  always  necessary  to  achieve  threshold  conditions  for  F? 
concentrations  of  0.3  percent  or  more.  In  contrast  to  the  situation  noted 
previously  for  XeF,  the  lack  of  lasing  at  higher  F^  concentrations  may  be 
attributed  to  absorption  at  the  laser  wavelength  due  to  F^-  It  should  be 
noted  that  the  threshold  measurements  for  the  high  F^  concentrations  also 
represent  an  upper  bound  due  to  absorption  of  the  stimulated  emission  by  F^ 
in  the  inactive  end  regions  of  the  discharge.  This  single  pass  loss  by 
absorption  is  approximately  10  percent  per  mole  percent  of 


SECTION  IV 


DISCUSSION 


The  measured  laser  threshold  power  densities  in  KrF  and  XeF  are 

considerably  smaller  than  those  reported  for  a shorter  active  discharge 
2 3 

length  ’ and  are  also  smaller  than  those  measured  in  the  present  device 
with  a short  pulse  (approximately  100  nsec)  discharge  supply.  Pump  power 
densities  for  several  different  XeF  experiments  are  plotted  in  Figure  5 
as  a function  of  £t/x,  where  t is  the  time  to  reach  threshold  measured  from 
the  start  of  the  discharge  pulse,  B,  is  the  active  gain  length,  and  x is  the 
cavity  round  trip  time.  The  pump  power  densities  are  the  sum  of  the  discharge 
and  e-beam  contributions.  The  e-beam  deposition  was  estimated  from  the 

4 

tables  of  Berger  and  Seltzer  and  then  multiplied  by  a factor  of  2.5  to 
account  for  multiple  scattering  effects.'* 


These  experiments  can  be  interpreted  as  follows.  At  1 atm  pressure 
the  primary  loss  mechanism  in  the  absence  of  stimulated  emission  is  radiative 
decay  for  lifetimes  of  the  order  of  10  nsec.  The  dominant  collision  partner 
for  deactivation  of  the  lasing  species  is  probably  F?  or  NF^,  and  collision 
cross  sections  in  excess  of  gas  kinetic  would  be  required  for  collisional  de- 
activation to  compete  with  radiative  decay.  Thus,  the  time  to  reach  thresh- 
old lasing  conditions  from  the  onset  of  the  discharge  excitation  is  determined 
by  the  number  of  round  trips  necessary  to  build  up  the  cavity  flux  so  that  it 
exceeds  the  radiative  loss  and  can  thus  begin  to  saturate  the  transition. 
Provided  that  the  mirror  losses  are  small  compared  to  the  round  trip  gain, 
the  delay  time  can  then  be  estimated  in  terms  of  the  small  signal  gain  from 
the  expression 


R exp[2(gQ  - a) «.] 


t/x 


where  0 is  the  half-angle  defining  the  internal  modes  that  are  stable  within 
the  cavity,  R is  the  net  round-trip  reflectivity  in  the  cavity,  including  all 
losses,  and  gQ  and  a are  the  small  signal  gain  and  absorption,  respectively. 
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□ Ref.  3,  R = 0.70 
v Ref.  6,  R = 0.97 


Figure  5.  Puinp  Power  Density  as  a Function  of  Delay  Time  for  XeF  Lasing. 

i he  continuous  curves  are  values  of  the  small  signal  gain  as  a 
function  of  £t/T  calculated  from  the  equation  presented  in  the 
text  for  9.  - 50  cm,  t = 8.7  nsec,  and  8 = 5 mrad.  The  scales 
for  the  small  signal  gain  and  pump  power  density  have  been 
arbitrarily  adjusted  with  respect  to  each  other  to  illustrate 
the  functional  dependence. 
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Since  the  small  signal  gain  is  directly  proportional  to  the  input  power 
density,  it  is  evident  that  threshold  is  defined  by  the  constant  product 
cf  the  input  power  density  and  the  delay  time  t.  Under  these  circumstances, 
threshold  conditions  are  effectively  determined  by  the  input  energy  density. 

To  test  this  analysis,  we  have  plotted  as  continuous  curves  in 
Figure  5,  values  of  the  small  signal  gain  gQ  as  a function  of  It/x  calcu- 
lated from  the  above  equation  for  three  reflectivities  and  a cavity  half 
angle  of  5 mil  1 i radians . The  scales  for  the  small  signal  gain  and  the  pump 
power  density  were  arbitrarily  adjusted  with  respect  to  each  other  to  illus- 
trate the  functional  dependence.  With  the  exception  of  the  electron-beam 
pumped  experiment  of  Brau  and  Ewing^,  the  data  follow  the  predicted  trend. 
The  implication  of  this  analysis  is  that  the  threshold  pump  power  densities 
measured  in  these  experiments  can  be  applied  to  other  experimental  systems 
using  the  above  relationships.  It  should  be  noted  that  for  a well  designed 
cavity,  increasing  the  path  length  does  not  change  Jl/x.  Thus,  the  threshold 
pump  power  is  independent  of  length  and  is  determined  by  the  maximum  per- 
missible delay  time  consistent  with  discharge  stability. 

Since  we  have  measured  the  delay  time,  t,  it  is  possible  to  estimate 
a small  signal  gain  in  the  present  experiments  using  the  above  equation. 

This  leads  to  small  signal  gains  of  the  order  of  1 percent/cm  and  thus 
satisfies  the  requirement  that  the  round-trip  gain  be  much  larger  than  the 
mirror  loss. 

It  is  also  possible  to  roughly  estimate  the  gain  from  the  input  pump 

★ 

power  density.  The  electric  discharge  energy  produces  metastable  A r with 

an  efficiency  of  approximately  25  percent  (Ref.  3).  About  75  percent  of  the 

* * 

Ar  reacts  with  Xe  to  form  Xe  , which  subsequently  reacts  with  NF,  to  form 
XeF  . If  radiative  loss  of  XeF  is  dominant,  then  at  steady  state  the  small 
signal  gain  is  approximately 
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g0  = o[XeF*]  = a 


1%/cm 


where  P = 15  kW/cm3  is  the  discharge  power  density,  n = (0.25)  x (0.75)  is 
the  product  of  the  efficiency  of  discharge  production  of  Ar*  with  the  branch- 


ing ratio  for  conversion  of  Ar*  to  Xe*,  = 11  eV/molecule  is  the  Ar*  ex- 
citation energy,  and  °T^ep*  = 7.5  x 10"24  cm2-sec  is  the  product  of  the  XeF* 
stimulated  emission  cross  section  and  radiative  lifetime  (this  product  is 


independent  of  rXep  since  o ~ VTXep*). 


In  the  case  of  e-beam  excitation,  the  energy  is  deposited  by  ioniza- 
tion of  the  Ar  diluent.  Foran  average  energy  loss  of  26  eV  per  ionization, 
the  estimated  e-beam  input  power  density  of  30  kW/cm3  corresponds  to  an  Ar 
production  rate  of  7 x 1021  Ar  cm'3-sec_1.  The  branching  ratio  to  form  XeF* 
is  not  yet  known;  however,  if  it  is  assumed  to  be  unity,  then  at  steady 
state,  the  small  signal  gain  would  be  approximately 


gQ  = o[XeF*]  - °[RTXeF*]  = (7.5  x 10“ 2 4 cm2 -sec)  (7  x 102 1 cm"3  sec-1 ) ~ 5%  cm"1 

This  is  much  higher  than  the  small  signal  gain  estimated  from  the  delay  time 
analysis,  which  implies  that  the  branching  ratio  for  Ar+  conversion  to  XeF* 
is  significantly  less  than  unity.  There  is  presently  not  enough  information 
about  the  steps  in  the  kinetic  chain  responsible  for  XeF*  formation  to  allow 
an  assessment  of  what  process  is  responsible  for  this  inefficiency. 

In  addition  to  the  above  threshold  power  density  measurements,  some 
experiments  were  made  with  30  percent  output  coupling  to  measure  the  pulse 
energy  and  efficiency  in  XeF.  The  experimental  conditions  for  two  measure- 
ments made  with  both  the  100  nsec  and  500  nsec  pulse-forming  network  power 
supplies  are  presented  in  Table  1.  Maximum  pulse  energies  of  approximately 
35  mJ  were  extracted  from  an  active  optical  volume  of  0.19  £ corresponding  to 
an  energy  density  of  0.18  J/£-atm.  The  electrical  efficiency  defined  by  the 
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Table  1 

XeF  Laser  Output  Energy  for  Two  Pump  Power  Densities 


Gas  Mixture 

Ar+1.2%  Xe+0.25%  NF3 

Ar+1 % Xe+0.1%  NF3 

Power  Supply 

100  NSEC  PFN 

500  NSEC  PFN 

Discharge  Voltage  (kV/cm-atm) 

3 

2 

Discharge  Power  (kW/cm3-atm) 

230 

60 

E-Beam  Power  (kW/cm3-atm) 

30 

30 

Discharge  Energy  (J/2.-atm) 

28 

20 

E-Beam  Energy  (J/Z-atw) 

5 

10 

Laser  Energy  (mJ) 

25 

35 

( J/H-atm) 

0.13 

0.18 

Electrical  Efficiency 

0.4% 

0.6% 
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ratio  of  lasing  energy  extracted  divided  by  total  energy  deposited  in  the  same 
volume  was  approximately  0.5  percent.  These  experiments  were  made  with  an 
applied  electric  field  of  approximately  2 kV  per  cm-atm  and  with  the  discharge 
input  power  density  a factor  of  2 larger  than  the  e-beam  power  density.  The 
lasing  pulse  lengths  of  these  experiments  were  approximately  100  to  200  nsec, 
the  discharge  becoming  unstable  after  approximately  300  nsec.  It  should  be 
observed  that  no  special  precautions  were  taken  to  optimize  the  discharge 
geometry  and  thereby  minimize  the  field  concentration,  as  might  be  achieved 
by  careful  shaping  of  the  discharge  electrodes. 

Based  on  the  present  results  for  the  threshold  pump  power  density,  it 
is  evident  that  efficient  laser  operation  in  KrF  and  XeF  could  be  achieved 
in  a 1-m  long  e-beam  sustained  discharge  with  input  power  densities  in  the 
range  of  50  to  100  kW/cm3  over  a 1-psec  time  scale.  Stable  discharges  have 
been  operated  in  KrF  and  XeF  at  these  power  densities  which  are  an  order  of 
magnitude  smaller  than  the  power  densities  used  in  short  pulse  (100  nsec  or 
less)  e-beam  sustained  discharges. 


« 
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